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Abstract. We prove an upper bound proportional to the surface area for the bipartite 
entanglement of the ground state and thermal states of harmonic oscillator systems with 
C*~) ■ disorder, as measured by the logarithmic negativity. Our assumptions are satisfied for some 

, standard models that are almost surely gapless in the thermodynamic limit. 

I. Introduction 

On the one hand, quantum computation and information processing depends in an essential 
way on entangled quantum states and, on the other hand, it is entanglement that is the most 
serious limiting factor in the numerical simulation of many-body quantum systems. Both 
Q-i! facts provide strong motivation to study entanglement of the ground states and equilibrium 
states of important systems such as oscillator lattice models. Furthermore, it was noted in [6] 
that entanglement implies a quantum contribution to black hole entropy. Consequently, the 
entanglement of ground states and temperature states of extended systems has been intensely 
studied, with establishing an area law bound as one of the main goals. 

Vidal and Werner [32] made the case for using the logarithmic negativity as a measure 
of entanglement and they showed how it can be computed in a number of cases, including 
Gaussian states for bosonic systems. The logarithmic negativity provides an upper bound 
for the widely used notions of entropy of entanglement in the case of pure states and the 
distillable entanglement in the case of mixed states. Following the ideas of Vidal and Werner, 
the logarithmic negativity has been calculated or estimated for a number of deterministic 
bosonic systems, primarily lattice systems of coupled harmonic oscillators [2J [2H CE21 US] • More 
recently, the logarithmic negativity has also been used to quantify entanglement in relativistic 
quantum field theories [TUj . 

In this paper we prove an upper area law bound for the logarithmic negativity for a class of 



cn 
o 



!h . disordered harmonic lattice models. The previous results of this type that are not restricted 



to one dimension, for either deterministic or disordered systems, all assumed the existence of 
a spectral gap above the ground state. For disordered systems, it is often no longer natural 
to suppose that there is a positive lower bound for the spectral gap uniform in the volume. 
Therefore, we will not make this assumption in this work. In the special case of one dimension 
an alternative approach may be possible following [7]. 

For this paper we have opted to include sufficient background information so that we can 
give complete proofs of results that often have appeared in a more restricted setting in the 
literature. E.g., in the Appendix we discuss the notion of partial transpose for operators on a 
the tensor product of two separable Hilbert spaces. We also included a proof of the result by 
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Vidal and Werner [32] that the logarithmic negativity is an upper bound for the entanglement 
entropy of a pure state in that setting. Similarly, in Section [3]the calculation of the logarithmic 
negativity of a quasi-free state is also discussed in reasonable generality. The main results are 
stated and explained in Section |2j 
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2. Harmonic Oscillator Lattices 

Our approach to proving entanglement area laws consists of two steps. The first reduces 
estimating an appropriate measure of entanglement for the oscillator lattice systems to specific 
properties of certain one-particle operators, known in the Anderson localization literature as 
eigenf unction correlators. We carry out the first step (Theorems 12.21 and 12. 3j) in a general 
context which we explain in Section 12.11 The second step requires a detailed analysis of the 
localization properties near the bottom of the spectrum of a one-particle system, which at 
present has been completed only in a number of specific examples, including the standard 
Anderson models. See Section 12.21 for a detailed description of these examples. 

2.1. Set-Up and Main Results. Consider a graph G = (I\£) where T is a countable set 
of vertices, often called sites, and £ is a set of undirected edges, i.e., pairs of vertices. We 
will assume that G is connected, i.e., for any x,y G T with x ^ y, there exists a path 
lx, y = (x = Xo,xi, • • • , x n — y) of finitely many vertices with (xj-i, Xj) G £ for all 1 < j < n. 
By d(x, y) we will denote the distance between x and y, defined as the minimum number of 
edges in a path connecting x and y and set d(x, x) = 0. In addition, we will assume that the 
graph is of bounded degree, i.e. 

(2.1) iV max :=sup|{?/Gr: (x,y) G £}\ < oo. 

This assumption is equivalent to the existence of a constant fi > such that 

(2.2) C„ = sup V e -^ d{x ' y) < oo. 

If the degree is bounded by iV max , it suffices to take fi > log A max for (12. 2p to hold. Conversely, 
if (12. 2p holds for some /i > 0, one has iV max < . 

For many interesting examples, G = £) for some integer v > 1 and edge set given by 
nearest neighbor pairs, e.g., with respect to the ^-metric, but this is just a special case. We 
also have results for graphs with exponentially growing volume such as the Bethe lattice. 

To formulate area laws, we need a notion of boundary. For this let Ao C A C T and assume 
that Ao is finite; A may be infinite. By <9Ao, we will denote the boundary of Ao which is given 
by 

(2.3) <9A = {x G A | there exists y G A \ A with (x, y) G £}. 

Although <9A depends on A, in general this dependence is of little importance and we suppress 
it in the notation. 

Given a graph G as above, we will consider oscillator systems defined as follows. To each 
site x G T, we associate a Hilbert space L 2 (M, dq x ) where we have used q x to denote the 
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spatial variable. It will be clear that our methods easily extend to the case where the single 
site Hilbert space is taken to be L 2 (R n , dq x ), and so we restrict our attention to the case of 
one- dimensional oscillators; mainly to ease the notation. For any finite set A C T, a Hilbert 
space "Ha is defined by setting 



(2.4) H A = L 2 (R, dq x ) = L (R , dq) 

with q = (q x ) x( z\. In each finite volume A and for any x G A, we will also use the notation 
q x to denote the position operator, i.e. the operator of multiplication by q x in "Ha, and by 
p x = —id/dq x we denote the corresponding momentum operator. By standard results, see 
e.g. p5], these operators are, on suitable domains, self-adjoint and satisfy the commutation 
relations 

(2.5) [q x , q y ] = [p x ,p y ] = and [q x ,p y ] = iS X)V l for all x, y G A . 

The models we consider will be defined in terms of two real- valued sequences {hx)y}x,yer 
and {^j/}x.yer- For any finite A C T, we will denote by the |A| x |A| matrix with entries 
P4 ?) kv = hit for all x,y G A and similarly Throughout this work, we will assume the 
following. 

Assumption 2.1. There exists a non-decreasing, exhaustive sequence of finite volumes A n C T, 
for n > 1, on which the matrices and nf 1 are real, symmetric, and positive definite. 
Moreover, we further assume that there exists C < oo for which 



(2-6) max \\hfj, \\{h®)-% \\h®\ 

uniformly in n. 



< C 



Note that when V itself is finite, this assumption may be applied to the constant sequence 
A n = T. We remark that, while requiring that the are positive definite and thus invertible, 

we do not assume a uniform bound on (/i^ In our applications this will amount to not 
requiring a robust ground state gap for the oscillator systems to be introduced next. 
Now, for any A n in a sequence satisfying Assumption 12.11 the formula 

(2.7) H n = Yl • i'-i'Zi',) 

x,y&A n 

= (q T ,p T )U n (q,p) 

defines a self-adjoint operator on H.A n which we refer to as a finite volume oscillator Hamil- 
tonian. In the final line above, we view q = (q x ) and p = (p x ) as column vectors indexed by 
x G A n with transposes q T and p T regarded as corresponding row vectors, and in this case, 
this line is a result of standard matrix multiplication with 



(2.8) H„ 







h A 



(p) 



When A n is understood to be fixed, we will just write and to ease notation. More 
general Hamiltonians could also be considered for the conditional statements made below. 
However, the only examples for which we can verify that the conditional statements hold are 
of the form described above, and so we will restrict our attention to this case. 
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It is well-known (e.g. [22] for more details) that H n can be written as a system of free 
Bosons, i.e., 

|An| 

(2.9) H n = ^2y e (2b* t b t + 1) , 

where the operators bi satisfy canonical commutation relations, i.e., 

(2.10) [b e , b e ,] = [b* e , b*,] = and [b e , b*,] = 5^1 for all £,£' G {1, 2, • ■ ■ , |A n |}, 

and the numbers 7^ > are the square roots of the eigenvalues of the positive definite matrix 

(2.11) K = {hto)V 2 hto\h®yl 2 . 

1 II 

For this reason, we often refer to h n as the effective single-particle Hamiltonian corresponding 
to H n , regarding it as a self-adjoint operator on £ 2 (A n ). 

In this work, we are interested in random models. In this case, we regard the components 
of sequences defining and as random variables on a probability space (f2,P). We will 
assume that there is a deterministic sequence of volumes A n for which Assumption 12.11 holds 
almost surely. By E(X), we will denote the expectation (average) of a random variable X on 
Q with respect to P. 

Our first result concerns the ground state of the random Hamiltonian H n . From the form 
of (12. 9p . it follows that, almost surely, H n has a unique, normalized ground state Q n G ~H\ n 
which is characterized by b(VL n = for all 1 < £ < |A n |. Let us denote by p n the orthogonal 
projection onto Q n . 

We study the bipartite entanglement in the ground state (and later also in the equilibrium 
states) with respect to a partition of the system. Explicitly, fix a finite set Ao C T and for 
any n > 1 large enough so that A C A n , write %h n = Hi <8> % where 

(2.12) Hx = (g) L 2 {R, dq x ) and H 2 = (g) L 2 (R,dq x ) 

zGAo x£A„\A 

Denote by p\ = Tr% 2 p n the reduction of the ground state projector to 7ii, and for any 
non-negative p with trace 1, let S(p) = — Trplnp be the von Neumann entropy of p. 

Theorem 2.2. Under the conditions stated above, in particular 112. fy) and Assumption \2.1\ 
and assuming further that there is a C < oo and a p! > p for which 

(2.13) E (K^) 1 / 2 ^,^ 1 / 2 ^) 1 ^)!) < C'e-M x >ti 
for all n > 1 and all x, y G A n . Then there exists C" < oo for which 

(2.14) E (S(pl l )) < C"\dA \ 
for all n > 1 . 

It has been argued that a bound on the entropy of entanglement of the type (I2.14p indicates 
that the ground state properties are computable. See the discussion in Section VI of [H] for 
an overview and [15] for a more nuanced discussion of the question in the case of oscillator 
lattices. 

Our proof of Theorem 12.21 uses the following bound for the entropy of the restriction of a 
pure state to one factor of a bipartite decomposition TL = Hi ® "H 2 : 

(2.15) S(Tr H2 p)<\og\\p T %=Af(p). 
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Here, the RHS is the log of the 1-norm of the partial transpose of the pure state p with respect 
to the given decomposition. This quantity is called the logarithmic negativity of the density 
matrix p with respect to the decomposition and is denoted by N(p). See the Appendix for a 
detailed discussion of partial transposes and, in particular, Lemma [A. 41 for a proof of (12.1 5ft . 

For any (3 > 0, the equilibrium state (aka thermal state) at inverse temperature /3 of the 
finite system with Hamiltonian H n is given by the density matrix 

As discussed in [32] the logarithmic negativity is a reasonable measure of the bipartite en- 
tanglement not only for pure states but also for mixed states (such as thermal states). In 
particular it is an entanglement monotone, see e.g. [32] for details. 

Theorem 2.3. Fix f3 > 0. Under the conditions stated above, assume further that there is a 
C < oo and a p' > p for which 

(2.17) E(|<(/i( p) ) 1 / 2 4,^ 1/2 tanh(/3(/i n ) 1 / 2 )(/i( p )) 1 / 2 ^)|) < C'e^'^ 
for all n > 1 and all x, y G A n . Then there exists C" < oo for which 

(2.18) E(M(p p , n ))<C"\dA \ 

for all n > 1, where pp t7l is the density matrix defined in 112.16]) and N(pp,„) is the logarithmic 
negativity with respect to the decomposition defined in 112.12]) . 

The proof of these results is given in Section |U which in turn relies on results we derive in 
Section [3] and the Appendix. 

2.2. Applications. Applications of Theorems 12.21 and 12.31 consist in verifying the localization 
bounds (I2.13P and ( I2.17P for the underlying effective single-particle Hamiltonian h n in concrete 
special cases of oscillator systems H n . 

Among available results in single-particle localization theory, bounds of the form (12.131) or 
( I2.17P correspond to strong forms of localization, which have only been rigorously established 
for the Anderson models (or models closely related to it). This leads us to consider the special 
case 

(2-19) H n = £ (±- P l + 9 -^A + £ A(g, - q y f 

xeA n ^ ^ (x,y)e£ n 

of (12. 7p . Here S n = {(x, y) G £ : x, y G A n }. The masses m, coupling parameters A and 
disorder parameter g are positive constants. 

Assumption 2.4. The {k x } x& r are independent, identically distributed random variables, which 
are absolutely continuous with bounded density p supported in [0, /c max ] for some /c max > 0. 

We thus have h^> — j^l and h^ q > becomes the Anderson model characterized by its qua- 
dratic form 

(2.20) (f,h®g)= Yl KM ~ f(x))(9(y) - 9(x)) + § £ kj(x)g(x), 

{x,y)e£ x& 

f,gE -^ 2 (r). Their restrictions to A n satisfy (I2.6P due to the boundedness of the k x and the 
fact that T is of bounded degree, in fact, 

(2.21) \\hf n || < 2AA max + -gk max . 
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Moreover, h n = -^/ijgj and the bounds ( T2~T3]) and (IXTTjl become equivalent to the existence 
of C < oo and fjf > fi such that, for all n > 1 and x,y £ A n , 

(2.22) E (|<^,/i- 1/2 5 y >]) < Ce"^") 
and 

(2.23) Ed^^^^tanh^^ 2 )^)!) < C'e^ 1 ^, 
respectively. 

Due to Assumption 12.41 it holds almost surely that k x > for all x G T. Thus, almost surely, 
the matrices h n are positive definite in any finite volume A n . However, as the support of the 
density p contains 0, there is no volume independent deterministic lower bound of the form 
h n > C > 0. A lower bound of this form is essentially what was used in [2U d2j to show that 
deterministic exponential decay bounds as in (I2.22p and (I2.23P hold, thus implying area laws 
as in Theorems 12.21 and 12.31 

In terms of the oscillator system Hamiltonian H n this means that an area law for the 
entanglement entropy of ground and thermal states was found to be a consequence of a robust 
ground state gap. This is seen by representing H n as a free Boson system (12. 9p . which shows 
that the ground state gap of H n is given by 

(2.24) 2min 7 ^ = 2mma(h 1 n /2 ) = 2 (min a(h n )) l/2 , 

see [22] for more details. 

A central goal of our work here is to show that in disordered oscillator systems it is not 
necessary for an area law to require a robust ground state gap, as long as averages are consid- 
ered on the left hand sides of (I2.22p and (12. 23 p . Using a term which was first proposed in a 
related context for disordered quantum spin systems in [TJ] (see also [IS]), we can argue that 
the localization properties of the single-particle operator h n lead to a mobility gap, which has 
consequences for the disordered many-body system H n similar to those of a robust ground 
state gap for a deterministic system. 

Single-particle localization bounds similar to (I2.22p and (I2.23P have been used in [22] to 
prove certain characteristics of many-body localization in disordered oscillator systems, such 
as dynamical localization in the form of zero-velocity Lieb- Robinson bounds, and exponential 
decay of ground state and thermal state correlations. Of particular interest in this context is 
the left hand side of (I2.22p . which, due to the absence of a robust ground state gap, presents 
an example of a singular eigenfunction correlator. Appendix A of [22J provides a detailed 
discussion of localization bounds for singular eigenfunction correlators, based on earlier results 
in the theory of Anderson localization (such as recently reviewed in |30j). In particular, this 
leads quite directly to the first of the following applications of our results. 

Theorem 2.5 (Lattice Systems). Let G = (Z",£) with edge set given by nearest neighbor 
pairs and A n = [— n, n] v PI Z 1 ', and let H n be the disordered oscillator system Ii2.19\) over G at 



fixed disorder g = 1, satisfying Assumption \2.4\ Then the ground state and thermal states of 



H n satisfy area laws of the form \2. 14\ ) and i\2.18\) . respectively. 

Proof. For G = 7L V the volumes \{y G G : \x — y\\ < n}\ grow polynomially in n. Thus (12. 2p 
holds for any /i > and, due to Theorems 12.21 and 12. 3[ it suffices to verify (I2.22p and (I2.23P 
for some fi' > 0. Both bounds follow as special cases from Proposition A. 3(c) in [22], as the 
functions tpi{t) = t~ 1 / 2 and tf2{t) = t~ x l 2 tanh(/3t 1 / 2 ) both have analytic extensions to the half 
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plane {z : Kez > 0} and for t G (0, oo) satisfy bounds of the form \<p(t)\ < Ct a for some 
a > —1. In fact, for ip 2 one can work with a = as tanh(/3t 1 / 2 ) ~ t 1 / 2 near 0. □ 

Without going into detail, we mention two natural directions in which Theorem 12.51 can be 
generalized: 

(a) We have considered = ell, but this can be generalized to larger classes of matrices 
h( p \ at least for the ground state case. Using ( 12. lip we see that the left hand side of (I2.13P 
has the form 

(2.25) n\^Ah {p Y\h^r l l\h^) l ' 4 5 y )\). 

One can prove exponential decay for this using the same arguments as in the proof of The- 
orem 12.51 as long as one has an a-priori exponential decay bound for the matrix elements 
of (/i^) 1 / 4 . This can be shown if the hP> are positive definite, diagonally dominant band 
matrices, satisfying the uniform norm bounds required in ( 12. 6p . In this case one can use the 
analyticity of tp(x) = x 1//4 in the right complex half plane to show exponential decay of matrix 
elements of (/z^) 1//4 , using arguments similar to those in, e.g., [TTJ. For example, one can 
choose hf^' = cl + 8T, where T is the next-neighbor hopping operator and 5 < c/(2v). 

Dealing with thermal states would require more work, as (/ i W) 1 / 2 ^ 1/2 tanh(/3(/ in ) 1 /2)(/ l (p))i/2 ! 
appearing in (I2.17|) . does not factorize into fractional powers of and hS q \ 



(b) It is also natural to ask if Theorem 12.51 extends to general graphs G = (T, S) as long as 
they have polynomially bounded volume growth, i.e. the sets A n (x) = {y : d(x,y) < n} grow 
polynomially in n (uniform in x). The crucial ingredient into the proof of Proposition A. 3(c) 
of [22], which we use above, is the well-known Lifshitz tail argument leading to localization of 
the single-particle Hamiltonians h n near E = 0. This requires to know a deterministic lower 
bound of the form E\ — E > C/n 2 for the ground state gap of the discrete Graph Laplacian 
on A n (x). We are not aware of a general result establishing such a bound for graphs with 
polynomial volume growth. But whenever it is known, area laws for the ground and thermal 
states will follow. 



In general, graphs of bounded degree have exponential volume growth, with the prototypical 
example given by the Bethe lattice. It is known that in this case localization proofs for the 
single-particle Hamiltonians h n require sufficiently large disorder g, see e.g. (TJ and [2]. In 
fact, the results of [2] establish that for the low-disorder Anderson model on the Bethe lattice 
the extended states regime may extend all the way to the spectral boundaries. However, if 
the disorder is sufficiently large, then we get area laws on general graphs of bounded degree: 

Theorem 2.6 (Large Disorder). Let H n be the disordered oscillator system A2.19\) on a general 
graph G = (T,£), satisfying h2. 1\) and Assumption \2.J\ If the disorder parameter g > is 



sufficiently large, then the ground state and thermal states of H n satisfy area laws of the form 
^2.1J^ and K2. 18\) . respectively. 

Proof. The main difference to the proof of Theorem 12.51 is that we now need to show (I2.22p 
and ( I2.23p . respectively, for some // > /i, with fi the constant from (12. 2p . As explained there 
this means that we need to show that we can choose /i' > logiV max . As a matter of fact, we 
will show that by increasing g one can choose // arbitrarily large. This will follow by well 
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established methods, but we will provide some detail as large disorder localization of singular 
eigenfunction correlators has not previously been discussed in the literature. 

For any subset A C T let h A be the restriction of the Anderson model (I2.20j) to A, and let 

(2.26) G A (x, y; z) = (S x , (h A - z)- l 5 y ) 

be its Green function. By adjusting arguments in Section 4 of [30] to the case of general 
graphs considered here one gets the following fractional moment bound (which is essentially 
already contained in pQ): 

For every s G (0, 1) there exists a constant C\ = Ci(s,p) such that 

(2.27) E(\G A (x,y;z)\ s ) < 9l e ~W '/<7iiWM*,i0 

9 s 

for all A C T, x, y G A, g > and zeC\l In fact, if A is finite, then one may choose z£C, 
including real values. 

Next we relate fractional moments to eigenfunction correlators: 

For every s G (0, 1) there exists C 2 = C 2 (s,p) such that 

/ \ / r \V(2-«) 

(2.28) E f sup \(S x ,u(h A ) X i(h A )6 y )\ 1 < C 2 [\g\ s jE(\G A (x,y; E)\ s dE 

for all finite A C T, x, y G A, g > and bounded open intervals Id. 

This can be proven following the ^-dependence of the arguments in Section 6 of [30], where 
the case g = 1 is considered. 

To conclude the proof of (I2.22p and (I2.23P we now apply these bounds to A = A n and, as 
before, write h n = h An , n = 1, 2, . . .. 

Choosing u(x) = f3~ 1 x~ 1 / 2 tanh(/3x 1 / 2 ) (which satisfies |u| < 1) and / = (0,E max ) : = 
(0, 2AA^ max + \gk nmy ) (which by ( 12 .2 1 f) almost surely contains the entire spectrum of h n ) and 
combining (T2T2Tj) and f l2~28|) yields 

(2.29) E(|(4,^- 1/2 tanh(/3^/ 2 )^)|) < (iC^C.E^)^ /c^) d ^y) ^ 

This proves (12 .23j) for g sufficiently large. 

Proving ( 12.221) needs a bit more work, as this requires handling a singular eigenfunction 
correlator. This can be done by the Riemann sum argument previously used in the proof of 
Proposition A. 3(b) in [22]. Decompose / = (0, i? max ) into 

(2.30) ^=(^>^y J = 1,2,..., 
and combine f[2T2Tj) and fl2T28j) to get 

(2.31) E(|(4,/^ 1/2 cyi) < C 2 Cp {jfj) 1/2 |^'l^j e-^ log(93/Cl7V — )d ^ y) . 

For s G (2/3, 1) the series in ( 12 .3 1 1) is summable, which concludes the proof of (12.221) . again 
for g sufficiently large. □ 

We conclude our discussion of applications by acknowledging that the types of disorder 
in oscillator systems which we have been able to handle is rather limited, essentially only 
covering Anderson-type (diagonal) randomness in h^ q \ Considering other types of disorder, 
such as random masses m or coupling constants A in (12.191) . is physically equally plausible, 
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but not enough is known about the localization properties of the associated single-particle 
Hamiltonians. Thus the applications provided here (as well as in [22]) motivate further studies 
of single-particle random Hamiltonians, with the goal of covering other phyically relevant cases. 

3. Gaussian States and Their Logarithmic Negativity 

The main goal of this section is to prove Theorem I3.4[ which provides a formula for the 
logarithmic negativity associated with an arbitrary finite volume A for a class of quasi-free 
states, including the ground and thermal states of the harmonic oscillator models introduced 
in Section [2j As was the case in all previous results of this kind [3 |2U [ITJ [12] , we start from 
the ideas in |32j . 

The logarithmic negativity is an upper bound for the entropy of entanglement. In the case 
of pure states (T = 0), the latter is the von Neumann entropy of the state restricted to the 
observables localized in A . The restriction of a quasi-free state is again a quasi-free state 
and this property makes it possible to essentially reduce the calculation to diagonalizing a 
one-particle operator. Calculating the logarithmic negativity means finding the one-norm of 
a partial transpose of the density matrix of the state or, equivalently, finding the norm of the 
partial transpose of the state regarded as a linear functional on the algebra of observables. 
The partial transpose of a quasi-free state, although in general not a state, is again a quasi- 
free (i.e. Gaussian) functional. This property makes is possible to find an explicit formula 
for the logarithmic negativity, see f !3.44p . In Section H] we prove the Area Law bound based 
on this formula, that is we prove Theorems 12.21 and 12.31 Since quasi-free functionals other 
than states, in particular, quasi-free functionals that are not positive, have not been widely 
studied, we provide in this section the necessary elements needed for the proof of Theorem 13.41 
in reasonable detail. 

The strategy of this section is as follows. Fix Ao C T finite. As is discussed in Section |2l 
for any finite A C T, both the ground and thermal states of the oscillator systems we consider 
can be expressed in terms of a density matrix p: 

(3.1) u(A) = Tr [pA] for any A £ B(H A ) ■ 

Here we have suppressed the dependence of the state u and the density matrix p on the finite 
volume A. For A with A C A, the logarithmic negativity of p is defined by 

(3-2) AT(p) = log(||p^||i) 

where p Tl is the partial transpose of p with respect to the decomposition as in ( 12.1 2p . The 
one-norm of p Tl equals the norm of the linear functional 

(3.3) w Tl (A) = Tr [p Tl A] for any A £ B(H A ) , 

which is well-defined exactly when ||p Tl 1^ < oo (proving the latter will be part of our argument 
below). Motivated by this relationship, we will start by studying the partial transpose of 
quasi- free functionals on the Weyl algebra (see Sections 13.11 and 13.21) . defined by 

(3.4) ou T ^W(f)) = e-^ I » 

in terms of a real, symmetric, positive definite matrix M. This functional need not be a state 
(it is not necessarily positive), but a version of Williamson's Theorem (see Proposition 13.21 
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below), implies that there exists a symplectic matrix S which diagonalizes M and therefore 

(3.5) co T ^W(Sf)) = J]V^I/WI 2 

3 

where Xj > are the symplectic eigenvalues of M; again, more on this can be found in 
Section 13.31 By explicit construction, we demonstrate in Section 13.41 the existence of a trace 
class operator p such that 

(3.6) u T ^W(Sf)) = TT[pW(f)} . 

Since there can only be one trace class operator satisfying this relationship (see Lemma 13.11 
below), we conclude that p is unitarily equivalent to p Tl . It is then straightforward to find 
the one-norm of p using its explicit form, and this yields the expression of the logarithmic 
negativity in Theorem 13.41 

3.1. On Weyl algebras and quasi-free functionals. We begin by introducing Weyl al- 
gebras, or CCR algebras, in the abstract setting. Here, we are brief and refer the interested 
reader to [8] for further background and details. Next, we describe quasi- free functionals on 
the Weyl algebra; these can be regarded as generalizations of the well-studied class of quasi- 
free states. The need for this generalization stems from our interest in logarithmic negativity. 
In particular, the partial transpose of a density matrix associated to a state on the Weyl 
algebra induces a functional that is not necessarily a state. 

Let T> be any real-linear space equipped with a non-degenerate, symplectic bilinear form a, 
i.e. a : T> x T> — > K. with the property that if a(f, g) = for all / G T>, then g = 0, and 

(3.7) a(f,g) = -a(g,f) for all f,g eT>. 

The Weyl operators over T> are introduced by associating non-zero elements W(f) to each 
/ ET> which satisfy 

(3.8) W(f)* = W(-f) for each / G V , 
and 

(3.9) W(f)W(g) = e- l ^ 2 W(f + g) for all /, g G V . 

As is proven e.g. in Theorem 5.2.8 [5], there is a unique, up to ^-isomorphism, C*-algebra 
generated by these Weyl operators with the property that W(0) = t, W(f) is unitary for all 
/ G V, and HW^,/) — 1|| = 2 for all / G D \ {0}. This algebra, commonly known as the Weyl 
algebra (also CCR algebra) over D will be denoted by W(P). 

Fix a real-linear space T> and the corresponding Weyl algebra W = W(T>). u is said to be 
a quasi-free functional on W if 

(3.10) = eW-bfrf) f or all / G P 

where r is a real-linear functional and s is a symmetric, real bilinear form on T>. 

It is clear that equation (I3.10p uniquely defines a linear functional on a dense subalgebra 
of W. Due to the form of (I3.10p . such functionals are also referred to as Gaussian. Not all 
these functionals are states, i.e. positive linear functionals on W, even if we assume they are 
continuous. In fact, it is well-known, see [31] . that a functional of the form (13.101) is a state if 
and only if 

(3.11) a(f,g) 2 <s(fj)s(g,g) for all /, g G V . 
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Observe that the inequality (13. lip above is equivalent to 

(3.12) 2a(f,g)<s(fJ)+s(g,g). 

3.2. Oscillator Model Examples. In this section, we briefly review the fact that the ground 
and thermal states of the harmonic oscillator lattices introduced in Section 12.11 are quasi-free, 
in the sense discussed above. We end this section with a discussion on corresponding partially 
transposed functionals. 

Let A C T be finite. The role of T> is played by the complex Hilbert space T>\ = £ 2 {A) with 
the symplectic form related to the inner product by 

(3.13) a(f,g) = lm[(f,g)} . 

The corresponding Weyl algebra Wa has a concrete realization: for each / G £ 2 {A), it is 
well-known that 



(3.14) W(f) = exp 



iJ2(MfU)]qj + Mf(j)]Pi) 

. jeA 

defines a unitary Weyl operator in B{TL A ) satisfying (13. 8p and (I3.9P above. Here, for each 
j G A, qj and Pj are the position and momentum operators introduced in Section 12.11 

The following basic fact is important for us. Here Bi(7i\) denotes the trace class operators 
on Ha. 

Lemma 3.1. If A G Bi(H\) and 

(3.15) Tr [A W{f)\ = for all f G £ 2 {A), 
then A = 0. 

Proof. This is a consequence of irreducibility of the Weyl algebra Wa in B(Ha) (e.g. Proposi- 
tion 5.2.4(3) of j8], the Fock space representation of the Weyl operators used there is equivalent 
to the representation (I3.14p when working in the Hermite function basis of "Ha = £ 2 (R A )) and 
von Neumann's double-commutant Theorem, showing that Wa is weakly dense in B(Ha)- 
The identity (I3.15P implies that Tr AC = for all C G Wa- This carries over to general 
C G B(Ha) due to the fact that C n ^> C implies Tr AC n — > Tr AC. Finally, use that 

(3.16) = sup \Tr AC\ 

C£B{H A ), ||C||<1 

to conclude that A = 0. □ 

Given an operator H^, as in (12.71) . denote by p\ the orthogonal projection onto the unique, 
normalized ground state of H^. A ground state functional oj\ on Wa is defined by setting 

(3.17) u A (W(f)) = Tr [ PA W(f)} for all / G £ 2 (A) . 

Here, and in what follows, we will regard the set A as fixed and simply write u and p. It will 
also be convenient to identify £ 2 (A) = £ 2 (A; C) with £ 2 (A; R) © £ 2 (A; R), i.e., 

(3.18) / G £ 2 (A; C) ~ f=( j^j ) G £ 2 (A; R) © £ 2 (A; R) . 
In this case, one calculates that 

(3.19) a(f,g) = (Jf,~g) where J= ( \ ~ 1 V 
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and (•, •) is the inner product on the direct sum. To ease the notation, we will just write f — f 
where this identification is to be understood. A well-known calculation, see e.g. [26], shows 
that 

(3.20) u(W(f)) = e-fam 

where M is the positive definite matrix 

(/^vV'W) 172 

o {hfr^hf{hfr^ 



(3.21) M 



and h\ = (h^ )^ 2 h^ ^{h^f 1 ) 1 / 2 is as in (12. lip . M, as above, is proportional to the real part of 
the ground state covariance matrix, i.e., the 2|A| x 2|A| matrix C given by 

(3.22) C = (cjk) where Cjk = Tr [pxjXk] and x = ( ^ ^ . 

In fact, one easily checks that 2C = M — iJ with J as in (I3.19p . 

Similarly, for any /3 > and A finite, a thermal state functional up on Wa is given by 

(3.23) = Tr[Pf>W(f)} = e~^ M ^ 

where pp is the thermal state density matrix, see e.g. (I2.16p . and the final equality above is 
again the result of a well-known calculation. Here 

{hf) 1 ' 2 coth(/3/ii /2 )^ 1/2 (^ ) ) 1/2 

(/i P) )~ 1/2 coth(/3^ /2 )/i A /2 (/i A P V 1 /2 



(3.24) M p 



and it satisfies 2Cp = Mp — iJ for the corresponding thermal state covariance matrix. It is 
clear that both u and up define quasi- free functionals on Wa in the sense of (13.101) . 

Using (13. lip , one also readily checks that the quasi-free functionals introduced above, i.e. 
both u and up, are states. In fact, let 



(3.25) R 



Ml' 2 
M 2 1/2 



where Mi and M 2 are, respectively, the upper left and lower right entries in the matrix M 
from ( 13.211) . It is clear that R is symplectic, i.e., 

(3.26) R T JR = J 

and therefore, R leaves the symplectic form invariant, i.e. 

(3.27) a(f, g) = (Jf, g) = (JRf, Rg) = a(Rf, Rg) 
In this case, 

(3.28) a(f,g) 2 = \(JRf,Rg)\ 2 < \\Rf\\ 2 \\Rg\\ 2 = (f,Mf)(g,Mg) 
and (EUD holds with s(/, /) = (/, Mf). With the relation 

2 

(3.29) coth(x) = ' 



e 2x - 1 ' 

it is clear that M < Mp, and thus the above argument proves that all of these functionals are 
states. 
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In fact, the quasi- free states for finite oscillator systems are always given by a density matrix 
on the Hilbert space as in ( 13. ip . 

As indicated previously, we are mainly interested in the logarithmic negativity associated 
to the above states. Motivated by calculations in Appendix |Aj see e.g. ( 1A.8I) . we define a 
partially transposed ground state functional by setting 

(3.30) v*(W(f)) = e~^ f) with M = [ \ \ ) M ( \ \ ) , 

M is as in (13.2ip . Ao C A fixed, and P is the diagonal matrix with 

-1 if x G A 
1 if x g" A . 

This functional is quasi-free, self-adjoint (oj Tl (W(f)*) = u Tl (W(f)) due to (I3.8P ). and nor- 
malized such that u Tl (l) = 1, but in general, it is not positive and therefore not a state. 
A partially transposed thermal state functional, cuj 1 , is analogously defined by replacing M 
above with Mp as in (13.241) . 

3.3. Diagonalizing quasi-free functionals. In this section, we return to the general setting 
of Section [37TI to discuss the diagonalization of quasi- free functionals. Recall that for any real- 
linear space T> equipped with a non-degenerate, symplectic form a, a quasi-free functional cu 
on the Weyl algebra W = W(T>) has the form 

(3.32) u(W(f)) = e fr (/)-3*(/.fl for aU / G V , 

where r is a real-linear functional and s is a symmetric, real bilinear form. As is the case in our 
examples, we will assume T> is finite dimensional. In fact, without loss of generality, we will 
assume T> = £ 2 (A) for some finite set A. By diagonalizing u, we mean finding an automorphism 
a of the Weyl algebra (or, equivalently, a unitary or anti-unitary transformation on Fock space) 
for which 

(3.33) u)(a(W(f))) = JJ e-^ Afel/(fe)|2 

k=l 

for some G R and n = |A|. This will enable us to calculate explicitly unitarily invariant 
quantities such as the entropy or the p-norm of the trace class operator associated with u. 

We start by noting that, without loss of generality, we may assume that r = 0. This is 
because for any such real-linear functional r, there exists g G T> such that r(f) = a(f,g). In 
this case, the automorphism a defined by 

(3.34) &(W(f)) = W(g)*W(f)W(g) = W(-g)W(f)W(g) = e~^W(f), 
where we have used the Weyl relations (13 .8p and (13. 9p . satisfies 

(3.35) io(a(W(f))) = e-^ fJ) . 

For the next step, we diagonalize the form s. This is done by choosing a basis in T> and 
representing s in terms of a 2n x 2n real, symmetric matrix M. A well-known version of 
Williamson's Theorem (see, e.g., [28]) then provides the existence of a symplectic matrix S 
which diagonalizes M; we state this as Proposition 13.21 below. As is proven e.g. in Theorem 
5.2.8 of [8], such an S induces an automorphism on W, in terms of which f !3.33|) is then clear. 



(3.31) 



P 
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Proposition 3.2. Let n > 1 and suppose that M is a real symmetric positive definite 2n x 2n 
matrix. Then, there exists a symplectic S such that 



(3.36) S T MS 



C 
£ 



where L is a diagonal matrix with entries A^ > for k = l,...,n. The numbers A& are 
the positive eigenvalues of iM 1 / 2 JM 1 / 2 ; they are also the positive imaginary parts of the 
eigenvalues of MJ. If, and only if, in addition, 

(3.37) M + iJ>0 
we have A& > 1 for all k = 1, . . . , n. 

Proof. The existence of a symplectic matrix S that diagonalizes M in the prescribed fashion 
is a special case of Wiliamson's Theorem [31], a particularly simply proof of which is given 
in [2E]- This includes the statement that the A& are positive. It follows from the proof of 
Theorem 3 in [28] that the A& are the positive eigenvalues of iM 1 ! 2 JM 1 / 2 . 
If the additional condition ( I3.37P is satisfied, we have 

(3.38) 0<S T (M + zJ)S= c) +U= (fl ~c) 

After reordering, the RHS is block diagonal with two-dimensional blocks of the form 

< 3 - 39 > (t 2 

which is non- negative definite if and only if A^ > 1. □ 

Remark 3.3. In the special case where M is block diagonal with two n x n blocks, i.e., of the 
form 

then the hermitian matrix iM 1 ! 2 JM 1 ! 2 takes the form 
(3-41) zMVVM^ =i l ° -Miyi 12 



mI /2 mI /2 



the square of which is 



(3.42) (zMV W / 2)2 = ( MpM&P 

Thus, in this case, the symplectic eigenvalues of M can also be found e.g. as the positive 

1/2 1/2 1/2 1/2 

square roots of the eigenvalues of M x M2M1 . Since M 2 ' MiM 2 ' and M 1 M 2} are both 

1/2 1/2 

similar to M 1 M 2 M X they have the same eigenvalues, and therefore, they could, as well, be 
used in determining these symplectic eigenvalues. 
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3.4. A formula for the logarithmic negativity. Now, we can combine the information 
of the previous subsections to prove the expression for the logarithmic negativity of ground 
and thermal states in Theorem 13.41 The essential observation is that M, see e.g (13 .30 p . is real 
symmetric and positive whenever M is, and therefore, the first part of Proposition 13. 21 applies. 
This symplectically diagonalizes M. The second part of the above-mentioned proposition does 
not apply, however, because 

, . / i o y o -i \ / o -i W i o 

[6A6) \o f J \ i o ) t \i o Jvop 

i.e., the partial transpose does not preserve the symplectic form J. We begin with a statement 
of the main result. 

Theorem 3.4. Fix A C T finite. For any finite A C T with A C A and H\ as in \2. 7| ] 

satisfying Assumption 2.1, we have that the logarithmic negativity associated to the ground 
state, respectively thermal state, and the decomposition in Ii2.12\) is given by 

(3.44) AT(p) = iTr[P + log(L- 1 )] 

where 

(3.45) L = mI /2 ¥M 2 ¥MI /2 , 

P + is the orthogonal projection onto the subspace where L < 1, P is the diagonal matrix with 

-l if xeA 
1 i/i6A\A , 

and for the ground state we have 

(3.47) M X = (^ ) ) 1/2 (^a)" 1/2 (^ ) ) 1/2 and M 2 = (hfy^h^ih^y^ , 
while for the thermal state at inverse temperature /3 we have 

(3.48) M x = (/ i i p) ) 1 / 2 (/ iA )" 1 / 2 coth(/3(/ iA ) 1 /2)(^) ) i/2 

(3.49) M 2 = (/ iA D) )- 1 / 2 (/ iA ) 1 / 2 coth(/3(/z A ) 1 /2)(^))-i/2 

with h K = (hf) i / 2 h [ «\hfyi 2 . 

We prepare the proof by a lemma in which we consider the Weyl algebra over the one- 
dimensional vector space D = C, denoted by W(C), and generated by 

(3.50) W(z) = exp (i[(Rez)q + (lmz)p\) = exp ( -^(za + za*) ), z E C, 



(3.46) P, 



XX 



V2 X 

where a = (q + ip)j \/2 and a* = (q — ip) / \/2. 

Lemma 3.5. Fix A > 0. There exists a unique, self-adjoint, trace-class operator p\ on L 2 ( 
for which 

(3.51) Tr [p x W(z)} = e "3 A l z l 2 for all z G C . 

Moreover, 

'1 if A > 1 
{ ifX<l. 



(3-52) \\p x \\i 



16 B. NACHTERGAELE, R. SIMS, AND G. STOLZ 

Proof. Denote by {|n) | n — 0, 1, . . .} the orthonormal basis of L 2 (R) given by the eigenvectors 
of the standard harmonic oscillator Hamiltonian a* a + |. In this basis a* and a correspond to 
creation and annihilation operators, 



(3.53) a*\n) = y/n + l\n + 1), for n > 0, and a\n) = \fn\n — 1), for n > 1. 

Using the commutation relation [a*, a] = 11, one readily verifies 

\z\ 2 ) m (n + m)\ 



(3.54) (n|W(z)|7i) = 



m>0 



2 m (m!) 2 n\ 



see e.g. (XII. 55) of [21] for more details. Note that the formal Taylor expansion used here 
can be justified by an analytic vector argument, e.g. Theorem 8.30 in [33]. Since A > 0, the 
number 

(3.55) a = a (A) = ^— ^, 

satisfies — 1 < a < 1 and therefore, the operator p\ defined by 

(3.56) p x \n) = (1 -a)a n \n) for n > 0, 

is clearly self-adjoint and trace class. A well-known calculation (e.g., see again [21]) using the 
identity 

(3 57) 1 _ ( n + m )l 

v ; (l-x) m+1 ^ m\n\ 

y ' n>0 



-X" 



gives 

(3.58) Tr[p x W(z)\ = (1 - a) ^ a n (n\W(z)\n) 



n>0 

, v J^^y- w (-|z| 2 ) m (n + m)! 
^ «J e 4 2^Z^ a 2 m (m!) 2 n! 

n>0 m>0 v y 



(-^I 2 ) m (n + m)\ . 

a 



a _ a ) e ^y 1-1*1 ; y 

/ -— ' 2 m (m!) ^ m!n! 

m>0 v y n>0 



e 4 



2\m 



2 m (m!) (1 - a) 7 



— g 4 g2(l-a) 

From the explicit form of p\, it is also clear that 

1 — a 

(3-59) UpaIK = - 



and therefore, we obtain (I3.52p . 

Finally, the uniqueness of p\ follows from Lemma 13.11 □ 



Proof of Theorem \3.4 



We verify f)3.44p for the ground state p. The argument for thermal states pp follows similarly. 
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According to L 2 (IR A ) = L 2 (R A °) ® L 2 (M A \ A °), Weyl operators W(f), f G f(A), can be 
decomposed as 

(3.60) W(f) = W(fV)®W(f®), 

where = /|a , = /|a\a - Thus, using flA.8j) . the partial transpose with respect to this 
decomposition becomes 

(3.61) W{f) T ^ = W(fU) T ® W(fW) = W ® W(fW) = W ( ( I M / 

where the last step uses the identification (I3.18P and the definition of P. It now follows from 
( E^DD that 

(3.62) Tr P W(fr = exp {-\ (( J ° ) /, M [ \ p ) /)) = 

where M is the real, symmetric, positive matrix defined in ( I3.30p . By Proposition 13.21 there 
exists a symplectic matrix S such that 



(3.63) S T MS 



£ 
C 



where £ is diagonal with entries Xj > 0, j — 1, . . . , |A|, the symplectic eigenvalues of M. With 
/ = Sg in (13.621) we get 

(3.64) TrpW^f 1 = e ~^ A ^ = Y[e~^ a 

j 

= HTr Px W( gj ) = Tr(^p Xj )W(g), 

j 

where p\. are the trace class operators on L 2 (IR) introduced in Lemma [3.51 

Letting U be the unitary that implements S in the representation of the Weyl algebra we 
are using (see [27] or [9]), i.e., 

(3.65) W(Sf) = UW(f)U* 
we observe that, for all g G £ 2 (A), 

(3.66) Tip W(Sg) n = Tr (® jPXj ) U*W(Sg)U = TrU (® jPXj ) U*W(Sg), 
or, as S is invert ible, 

(3.67) TrpW(f) Tl = TrU (® jPXj ) U* W(f) for all / € £ 2 (A). 

Now Corollary IA.2I proves that p Tl = U{®jp Xj )U*. In particular, p Tl is trace class and 

(3-68) l|p Tl ||i = IIlK-lli- 

j 

Using (I3.52p to calculate the trace norm of p Xj we get 

(3-69) A/-(p) = log H^H^^logllpAjl^ l °Zy- 

j y.\ r : j 
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By Remark 13.31 (applied to M) the Xj are the positive square roots of the eigenvalues of L 

1 

2 



Ml l2 ¥M 2 ¥Ml /2 . This shows that the right hand side of fl3~69l) is equal to \Ti[P + log(L -1 )]. 



□ 

4. An Area Law for the Logarthimic Negativity 

The goal of this section is to prove Theorems 12.21 and 12.31 This is done after establishing 
two deterministic facts. First, in Lemma 14.11 we prove an upper bound on the logarithmic 
negativity associated to ground and thermal states of oscillator systems. The arguments here 
follow [12] rather closely, however, we avoid making assumptions on the spectral gap of the 
one-particle operators involved. Next, we prove a simple geometric fact about the graphs we 
are considering in Lemma 14.21 The proofs of our main results follow. 



Lemma 4.1. Fix A C T finite. For any finite A C T with A C A and if a a s in (2.1) 



satisfying Assumption \2.1\ we have the following bound on the logarithmic negativity associated 
to the ground state, respectively thermal state, and the decomposition in i\2.12\) : 

(4.1) >V(p)<2||J^r 1 ||^ |(^,M 2 "\)|, 

xeA y£A\A 

where for the ground state we have set 

(4.2) Mf 1 = (/4 p y i/2 (/i A ) V2 (/4 P) r 1/2 and MT 1 = {hf) 1 ' 2 ^)- 1 ' 2 ^) 1 ' 2 
while for the thermal state at inverse temperature (3, 

(4.3) Mf 1 = (/ i i p y i / 2 (/ iA ) 1 / 2 tanh(/3(/z A ) 1 / 2 )(^ ) )- 1 / 2 

(4.4) M 2 -! = (^ ) ) 1 / 2 (/ iA )- 1 / 2 tanh(/3(/z A ) 1 /2)(^) ) i/2 

with h A = {hfy^hfihf) 1 ' 2 in both cases. 

Proof. Recall the results of Theorem 13 A[ in particular, (I3.44p and (13 .45 p . Using the concavity 
of the logarithm, we immediately obtain the following bound: 

(4.5) 2Af(p) = TrP+logL" 1 < TrP + (L" 1 - 1) . 

Next, we note that L~ l = Af~ 1/2 PAf 2 " 1 PM~ 1/2 can be written as L~ l = A + B with 

(4.6) A = M'^M^M^ 2 and B = M^PfMf 1 , P]M~ 1/2 . 

We claim A < 1. In the ground state case, we have M 2 = Mf . This directly implies A = 1, 
and there is nothing to prove. In the thermal case, recall that Mi and M 2 are given by (13.481) 
and (I3.49p . respectively. As a result, they are both positive matrices; hence, so too are L -1 
and A. It is clear that 

(4.7) M^M^ 1 = (/i A 0) )- 1/2 tanh 2 (/3(/ iA ) 1 / 2 )(/i A 0) ) 1 / 2 
and thus A is similar to tanh 2 ((3h 1 / 2 ) } i.e., 

(4.8) A = i^tanh 2 (/3/i A /2 )J^- 1 with K = M\' 2 {hfy 1 ' 2 . 

We conclude that the spectrum of the positive operator A is contained in [0, 1], and hence we 
have A < 1 as claimed. Using this fact we obtain 

(4.9) 2Af(p) < TiP + (A- t + B) < TiP + B. 
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To proceed, note that 

(4.10) Tr P + B = Tr P + M[ 1/2 ¥[M^ , P] M~ 1/2 P + 

< \\P + M^ 1/2 \\ 2 WFlM^^Wi 

< \\p + Mr 1/2 \\ 2 E k<^p[m 2 -\p]^)|. 

x,ydA 

The first factor is of course bounded by 1 1 ikTj - 1 1 1 . For the second factor we used the fact that 
the 1-norm can be bounded by the sum of the absolute values of the matrix elements in any 
basis. Observing that 

/a 1 1 \ /r p r. r i pi A \ / if x,y G A or x,y G A\ A , 

(4.11) (5 X , P [M 2 , Pj 5 y ) = | _ 2(4; M _ l6y) otherwise . 

leads to the claimed result. □ 

As observed in [12], if h\ > c > 0, P + = for sufficiently small /3, meaning that all en- 
tanglement vanishes at high temperatures. The previous result holds, however, at all positive 
temperatures. 

For the sake of completeness, we prove a general fact about the graphs we are considering. 
Lemma 4.2. Let A C T be finite and suppose fi> is such that 



(4.12) LV=sup^e 



-fj,d(x,y) < ^ _ 



Then for any A C A C T, 

(4.13) E ^ d{x ' v) < Cl\dA \ 

xGAq yeA\A 

where \dA Q \ is as defined in Ii2.3\) . 

Proof. Note that for every x G A and y G A \ A there is at least one u G <9A such that 
d(x, u) + d(u, y) = d(x, y). Therefore 



(4-i4) E E ^ d[x>y) ^ E E 



e -fid(x,u) e -fid(u,y) 



X&A y£A\A uGdA x£A ,y£A\A 

d(x ,u)-\-d(u,y) — d(x ,y) 



u&8A x£T,y£F 



<C^\dA \. 



□ 



We can now complete the proofs of Theorems 12.21 and 12.31 
Proof of Theorem \2.3[ Using Lemma 14.11 and Assumption 12.11 it is clear that 
(4.15) E(M( P ))<CJ2 E E(|(^,M 2 -\}|) 

x£A y£A\A 

The quantity on the right-hand-side above is precisely what appears in ( 12. 17ft . ( 12.181) now 
follows from Lemma [4.21 □ 
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Proof of Theorem \2.2 . The proof of Theorem 12.21 follows as above, after an application of 
Lemma IA.4I □ 

Appendix A. On Partial Transposes and an Entanglement Bound 

The goal of this section is to collect some basic facts about partial transposes valid in 
the context of infinite dimensional complex Hilbert spaces. We begin with a discussion of 
general conjugations and their corresponding transposes in Section IA.ll An example of the 
transpose of a Weyl operator with respect to the natural conjugation ends this section. Partial 
transposes are defined in Section IA.21 and a number of important properties, each used in 
the main text, are proven in detail. Finally, Section TA.3I contains a proof, originally given by 
Vidal and Werner in [32], that the logarithmic negativity provides an upper bound on the 
von-Neuman entropy of the restriction of a pure state. 

A.l. Conjugations and Transposes. First, recall the definition of the operator transpose 
with respect to a given conjugation in a Hilbert space. 

Let (H, (•, •)) be a separable complex Hilbert space and C : Ti — > H be a conjugation, i.e. 
(Cf, Cg) = (g, f) for all f,g£H and C 2 = I. 

This implies that C is anti-linear, i.e. C(af + (3g) = aCf + (3Cg for all a, j3 G C, /, g G H. 

An important fact is that C is a conjugation in % if and only if there exists an orthonormal 
basis (ONB) {e^} of "H such that 



\ k / k 



(A.l) C \^ x k e k = x k e k for all (x k ) G £ 2 



To see this, first note that for an anti-linear mapping C the property flA.ip is equivalent to 
the existence of an ONB such that 

(A. 2) Ce k = e k for all k. 

It's easy to see that (IA.20 is sufficient for C to be a conjugation. On the other hand, if C 
is a conjugation, then by Zorn's lemma there is a maximal orthonormal system {e k } with the 
property flA.2j) . To show that this is an ONB consider the orthogonal complement D of the 
subspace spanned by {e k }. If D were non-trivial, then it would contain a normalized vector u. 
One can now construct a normalized vector v G D with Cv = v, contradicting the maximality 
of {e k }. This is done separately for the two cases Cu G span{u} and Cu G" span{w}. In the 
first case v can be chosen as a suitable scalar multiple of u, in the second case one may choose 
v = (u + Cu)/\\u + Cu\\. 

While many different conjugations exist on any given Hilbert space, it follows that any two 
conjugations C and C on % are unitarily equivalent. In fact, for an ONB {e k } such that 
Ce k = e k , let U be the unique unitary such that 

(A.3) Ue k = e k 

for all k. Then it is easily seen that 

(A.4) C = U*CU. 

For any A G B{T-L), the bounded linear operators on H, its transpose A T with respect to 
the conjugation C is defined by 

(A.5) A T = CA*C. 
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We have A T £ B(H) with ||A T || = ||A||. Other basic properties of transposes which we will 
use below are (A*) T = (A T )* and that U is unitary if and only if U T is unitary. 

With an ONB {e^} associated with C via ( 1A.1I) (and only with such ONBs), the transpose 
with respect to C is characterized by 

(A. 6) (e,-, A T e k ) = (e k , Aej) for all j, k. 

For two conjugations C and C which are unitarily equivalent via U as in flA.4j) . a calculation 
shows that the corresponding transposes are unitarily equivalent via U T U, 

(A.7) A f = (U T UyA T (U T U). 

Example. If Ti — L 2 (X,fi) for a measure space (X, /i), then the natural conjugation on Ji is 
given by Cf = f and the associated transpose by A T f = A*f. An ONB {e^} of L 2 (X, /i) 
satisfies ( 1A.ll) if and only if all e k are ^-almost everywhere real-valued. In the special case 
Ti = ^ 2 (r) for a countable set T, this is true for the canonical basis e 3 -, j £ T. 

If A is finite and U = Ha = L 2 (R A ,dq) as in ([231) , then for each / e ^ 2 ( A ) the We y! 
operators are defined by ( 13. 14ft . The Weyl operators are unitary and their transposes with 
respect to natural conjugation are given by 

(A.8) (W(f)) T = W(f). 

This is seen by combining the Weyl relations ( 13. 9 p with the facts 

(A.9) 0(gf = 0(g), 0(pf = 0(-p), 

where, as usual, (p(q) and 4>{p) denote multiplication operators in q and p space, respectively. 
In particular, <fi(p) = F~ 1 (j)(q)F for the Fourier transform (Ff)(y) = (2tt)~ 1 ^ 2 J e~ lxy f(x) dx. 
Verifying the second one of the identities (1A.9[) uses that FC = RCF, where Rf(x) = f(—x). 



A. 2. Partial Transposes. We will also need the concept of a partial transpose of suitable 
classes of linear operators on tensor product spaces. For this, let "Hi and "H 2 be separable 
complex Hilbert spaces and W = Hi ® H-i- Let C be a conjugation on Hi and A H- A T 
transposition with respect to C on B(Hi). 

For tensor products A <g) B with A £ B(Hi) and B £ B(H 2 ), we define a partial transpose 
with respect to the first component of the tensor product as 

(A.10) (A® B) Tl := A T ® B = CA*C ® B. 

For finite dimensional Hilbert spaces Tii and %2 on e can uniquely extend the definition 
of partial transpose to general S £ B(H) by imposing that the mapping S >->■ S Tl is linear 
from B{H) to B(H). However, in the cases relevant here the Hilbert spaces have infinite 
dimension. In this case, an attempt to linearly extend (1A.10|) to general S £ B(9i) will have 



to allow unbounded partial transposes 5* Tl . This is due to the fact (e.g. jl]) that in the finite- 
dimensional case the norm of the mapping S t— > S Tl as a linear operator on B{T-L) is given 
by min(dim'%i, dim I-L2) and thus increases with the dimensions. From this it is not hard to 
construct an example of a bounded operator on an infinite dimensional space with unbounded 
partial transpose. 

We can avoid dealing with these issues here as in all our applications partial transposes will 
only have to be considered for operators which fall in one of two special classes of bounded 
operators, whose partial transposes are easily seen to be bounded. One of these classes are 
products A <g> B, which are covered by (1A.10B . so that \\(A <g> B) Tl \\ = \\A T <g> B\\ = \\A ® B\\. 



22 B. NACHTERGAELE, R. SIMS, AND G. STOLZ 

Another convenient class are the Hilbert-Schmidt operators on Ji, which we will denote by 
B 2 (H). Thus, let S G B 2 (H) and {e^} be the ONB of U x associated with C via fTATTj) . Also, 
let {fk} be any ONB of T-L 2 . Define a linear operator S Tl which acts on the basis vectors 
ej <g> fk of K as 

(A.ll) S Tl ej ®f k = ^(ej ® f t , Se m <g> / fc )e m ® 

Note that the meaning of partial transpose is reflected in the fact that its matrix elements are 

(A. 12) (e m ® f e , S Tl ej ® f k ) = (e,- ® f t , Se m ® 

which should be compared with (1A.6I) . Up to rearrangement, S Tl has the same matrix-elements 
as S, meaning that the Hilbert-Schmidt norm is preserved under partial transposition: 

(A.13) W^Wl = Y,\\S Tl e^f k \\ 2 = \(^®fe,Se m xf k )\ 2 

j,k j,k,m,£ 

= ^||^e m ®M| 2 = \\S\\l 

m,k 

In particular, the operator S Tl defined by flA.lip on a tensor product basis has a unique 
extension to a bounded operator on %. 

One can check that the definitions (1A.10I) and (1A.11I) are consistent, meaning that for 
S = A ® B, A G B 2 ('Hi), B G -82(^2), both definitions of S Tl coincide. Another property 
which holds for both classes is that (S Tl ) Tl = S. 

For the next result we denote the trace class operators on % by Bi(H). 

Lemma A.l. (a) Let S G B X (H) such that S Tl G B^U). Also, let A G B(Ui) and B G 
B(H 2 ). Then 

(A. 14) TiS Tl (A(g) B) = TrS(A® B) Tl . 

(b)IfS,Re B X (H) are such that 
(A. 15) Tr S (A ® B) Tl = Tr R (A ® B) 

for all A G B(Hi) and B G B(^-L 2 ), then R = S Tl and, in particular, S Tl G B 1 ('H). 

Proof, (a) Expanding Aei and Bf m with respect to the bases (e 3 -) and (/&), respectively, we 
get 

(A. 16) (e t ® / m , S Tl (A g) £?)e £ <g> / m ) = ^(e„ Ae £ ) (/ fc , Bf m ) (e e ® / m , S^e,- ® / fc ). 

i.fc 

Using flA. 12[) as well as (e^Ae^) = (ej, (A T )*ee), we see that the right hand side of ( 1A.16I) is 
equal to 

(A.17) (((A T )*e e ) <g> / m , Se £ <g> (S/ m )) = (e £ ® f m , (A T <g> /)£(! ® 5)e £ ® / m >. 

This allows to conclude by using cyclicity of the trace, 

(A.18) Tr S Tl (A®B) = Tr (A T ® I)S(I ® B) = Tr S(A T ® B). 

(b) Note that it is not a priori clear that S Tl G Bi(H). However, as S G Bi(7i) C B 2 (H), 
we know from the above that S Tl G B 2 {T-L). Thus it is a bounded operator on H. and therefore 
characterized by its matrix elements (e k ® f m , S Tl ej <g> ft). 
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Choosing A = \ej)(ek\ and B = \fe)(f m \ we find that Tr R(A ® B) — (e k tg) f m ,Rej <S> ft), 
while Tr S{A ® = (e fe <g> / m , S^e,- <g> //). We conclude from flATS]) that R and 5 Tl have 
the same matrix elements and thus coincide. □ 

For the special case of L 2 -spaces, it suffices to verify flA.15j) for Weyl operators: 

Corollary A. 2. LetUi = L 2 (R Al ) ; U 2 = L 2 (M A2 ) ; H = L 2 (R AlUA2 ) ; and let C ^ C Tl denote 
the partial transpose with respect to natural conjugation in the first component Hi ofH. If 
S, R G B\ ("H) are such that 

(A. 19) Tr SW(f) Tl = Tr RW(f) 

for all f G £ 2 (Ai U A 2 ), then R = S Tl . 

Proof. This can be reduced to Lemma [A. 1( b) using the arguments in the proof of Lemma [37T1 
First, (1X191 implies TrS(Wi ® W 2 ) Tl = Tr #(Wi <g> W 2 ) for all Wi G W Al , W 2 G Wa 2 . 

As W A . is weakly dense in B^H^f), j = 1,2, for all ^4 G Btf-L^A, B G BiTi^A there are 

W^- )T i G Waj, j = 1, 2, n = 1, 2, . . ., such that Wi |Tl — > A and H^ 2 ,n B. Weak convergence is 
conserved under transposition and tensor products, i.e. 

(A.20) W ltn (8) W 2 ,n ^ A <g> 5, W\ T n (8) ^> A T <g> 5. 

We conclude that ( 1A.15I) holds. Thus Lemma [A. 1( b) completes the proof. □ 

Unitary equivalence of transposes (1A.7I) with respect to different conjugations extends to 
partial transposes: 

Lemma A. 3. Let C and C be conjugations in Hi, which are unitarily equivalent via U as 
in (A. 4 )- Denote the corresponding transpose maps in B{U X ) by A^ A T and A ^ A T . Let 
S G B 2 (H) or S = A®B , A G B(Hi), B G B(H 2 ). Then S Tl and S Tl are unitarily equivalent 
via V = (U T U) <g> I, 

(A.21) S fl = V*S Tl V. 

Proof. It is very likely that flA.21j) will hold for larger classes of operators S if one goes through 
the effort of properly extending the definition of partial transpose, but it is most convenient 
here to merely check (1A.21I) separately for the two cases covered by flA. 101) and flA.llj) . 

For tensor products S = A <g> B, (IA.21I) follows immediately from (IA.7I) . 

For S G B 2 {H) we need to verify that 

(A.22) (e m <g> h, S% <g> f k ) = (e m <g> f £ , V*S n Ve j <g> f k ) 

for all m, t, j, k, using the ONBs introduced above. By (1A.3[) and (1A.12j) the left hand side 
of (1A.22|) is equal to 

(A.23) (ej ® f e , (U ® I)S(U* ® J)e m ® /*). 

For the right hand side of flA.22j) one gets 

(A.24) {{U T e m )®ft,S T ^{U T e j )®f k ) = ^ (e m , Ue r )(e 3 , Ue s )(e r <g> S Tl e s <g> / fc ) 

= 2^(e r , U*e m )(e s , U*ej)(e s <g> jfr, Se r g) 
= ((^ , e i )®//,5 , (CT , e TO )®/fc)- 
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The latter coincides with ()A.23j) . 



□ 

A consequence of Lemma IA.3I which we will use later is that S Tl G Bi(H) if and only if 
S Tl G Bi(H) and that in this case 

(A.25) \\S Tl \\i = ll^lli. 

A. 3. Entanglement bounds for pure states. In this final section, we prove - mainly for 
the sake of completeness - a result from [32], see Lemma IA.4I below. 

Let us briefly recall some basic definitions. Let Hi and H2 be separable complex Hilbert 
spaces and set H = Hi <S> H2- For any normalized vector Q G H, denote by p the orthogonal 
projection onto Q. By pi, we will denote the restriction of p to Hi, i.e. 

(A.26) Pl = Ti n2 p 

In the result we prove below, two quantities are of interest. The first, S(pi), is the von 
Neumann entropy of this restriction: 

(A.27) 5(p 1 ) = -Trp 1 log( Pl ). 

The second is the logarithmic negativity of p with respect to the above decomposition. To 
define it, we chose a conjugation on Hi and set 

(A.28) A/-(p) = log(||p T 1i) 

where p Tl is the partial transpose of p with respect to this conjugation. By flA.25j) this 



definition is independent of the choice of conjugation. Note that if p Tl ^ Bi(H), then the 
above is interpreted as infinite. The main result is as follows: 

Lemma A. 4 (Vidal- Werner) . Under the assumptions above, 

(A.29) S(pi) < Af(p) . 

Proof. We begin by writing Q in terms of its Schmidt decomposition, i.e., 

(A.30) ft = ^c Q |e a ®/ Q > 

a 

where each of {e a } C Hi and {//?} C H2 are orthonormal sets indexed by a common, countable 
set. Without loss of generality, we will take c a > 0. By normalization, it is clear that 
Yla c a = 1- Both quantities of interest in flA.291) can be calculated in terms of these coefficients 
c a . 

The restriction pi can be computed as 
(A.31) pi = y^ y c 2 a \e a )(e a \ 

a 

from which it is clear that the c 2 a are the eigenvalues of pi. In this case, 
(A.32) S(Pi) = -E c X c ') = 2 E ln f^) c ' 

a a x ' 
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the inequality following from Jensen's inequality. Without further assumptions, it is possible 
that the right hand side is infinite. 

For the logarithmic negativity, we only consider the case that ||p Tl ||i < oo; otherwise the 
bound ( 1A.29I) is trivial. Extend the orthonormal set {e a } C "Hi to an orthonormal basis, select 
a conjugation C leaving this orthonormal basis invariant, and denote by T\ the corresponding 
partial transpose. As remarked above, we may use any partial transpose to calculate the 
logarithmic negativity. We will use T\. 

As is clear from (1A.12[) . we have that 

(A. 33) p Tl = ^2 c aCp\ep <g> f a )(e a <g> f p \ . 

Now, introduce a linear mapping F : Ti — > % by declaring that 
(A.34) F(\e a ®fp)) = \ep®f a ) 

for all a, (3 corresponding to vectors in the Schmidt decomposition and extend by the identity 
otherwise. Then F is a unitary and furthermore, 

(A.35) p Tl = ^c a cp\ep® f a )(e a ® fp\ = ^ c a c p F (|e Q g> f p )) (e a ® f p \ = F(C X ® C 2 ) , 

a,l3 a,0 

where we have set 

(A. 36) Ci = ^2c a \e a )(e a \ and C 2 = ^ cp\fp){fp\ . 

a p 

We have then that 



(A.37) 



0^)11^ ||d ^C^HICillx-ll^l^ \J2< 



The bound claimed in (1A.29I) now follows from (IA.32]) . 



□ 
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